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We show that classical-quantum correspondence of center of mass motion in two coupled delta- 
kicked rotors can be obtained from intrinsic decoherence of the system itself which occurs due to 
the entanglement of the center of mass motion to the internal degree of freedom without coupling 
to external environment. 
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Classical-quantum correspondence in a classically 
chaotic system has been one of the most interesting prob- 
lem in physics for a long time In quantum me- 
chanics, the time evolution of a wave function follows a 
linear Schrodinger equation, and so there is no possibil- 
ity of a sensitivity on the initial condition, a trademark 
of classical chaos. Also, chaotic diffusion is suppressed 
by quantum localization ||. It has been revealed that 
some crossover time t r = log(I/h)/X (I is a characteris- 
tic action and A is a Lyapunov exponent) exists so that 
classical-quantum correspondence breaks down for t > t r 
j| . As h — ► the crossover time increases indefinitely 
and classicality is recovered. However, the problem is 
that h is a nonzero constant and t r is not sufficiently 
large considering its logarithmic dependence of h [||. 

Recently, the relation between decoherence and the 
classical-quantum correspondence has been investigated 
extensively Pj3|-|l2|] . Decoherence breaks the purity of 
initial superposition, which should be conserved in the 
absence of coupling to the environment, and thus only the 
partial fraction of whole Hilbert space, namely pointer 
states, are selected by the environment @]. The dynam- 
ics of the system coupled to the environment shows the 
unique characteristics of the system independent of the 
coupling strength as long as it is not too large or small 
p0| . In other words, with appropriate coupling to envi- 
ronment, the Lyapunov exponent or entropy production 
rates, which are important physical quantities charac- 
terizing a chaotic system, can be reproduced quantum 
mechanically. However, there has been some debate on 
whether decoherence from environment is indispensable 
to obtain the classical-quantum correspondence for clas- 
sically chaotic systems ||. 

In this letter, we show that decoherence can occur nat- 
urally in composed system even when we ignore the cou- 
pling to the outer environment. The initially pure center 
of mass states become dynamically entangled to the in- 
ternal degrees of freedom, which effectively acts like the 
environment. Let us consider a classical object governed 
by a Hamiltonian H = P 2 /2M + V(X), where M is the 
mass of the classical object. Since classical objects are 
composed of many particles, a complete Hamiltonian will 



be given by H = J2ipf/ 2m i + v x{ x i)] + Vi {x x , x 2 , ■ ■ •), 
where pi and Xi are momenta and coordinates of the 
constituents, respectively. We assume that the mass 
of each constituent m, is equal to m. If the force 
fi = —dVi(xi)/dXi is linear, we can ignore the motion of 
the individual constituent particles to describe the cen- 
ter of mass dynamics of the macroscopic classical object 
since MX=J2f l = kJ2%i = NkX, where X, N, and 
k are a position of center of mass, the number of con- 
stituents, and a constant characterizing the linear force, 
respectively. If the force fi is nonlinear, however, the 
description of the motion of the classical object by us- 
ing only Hq is not immediately clear. We already know 
from our experience, nevertheless, that the classical dy- 
namics of the macroscopic object is well described by the 
Hamiltonian Hq. It means that ^2 fi should almost co- 
incide with F(X) — —dV(X)/dX in values. Thus, we 
can ignore the difference between H and Hq to consider 
the motion of classical objects. However, the minor dif- 
ferences between H and Hq play an important role in 
quantum mechanics since the internal degree of freedom 
can induce decoherence on the center of mass motion. In 
this case, the correct classical-quantum correspondence 
of the center of mass motion will be obtained not from 
Hq but from H. 

Only due to computational difficulty, in this letter we 
consider a two particle system. Even with this sim- 
ple model, we have observed some evidence of the in- 
trinsic decoherence originating from the entanglement to 
internal dynamics. We will show the delocalization of 
wavefunctions in momentum space and the coincidence 
of classical and quantum entropy production rates. In 
the previous study of two coupled quantum kicked tops, 
the entanglement of two initially decoupled kicked tops 
increases linearly in time with a rate which is a linear 
function of the sum of the positive Lyapunov exponents 
However, the classical limit has not been taken, and 
thus the classical-quantum correspondence has not been 
considered although signature of classical chaos has been 
observed for low quantum numbers. A model of two in- 
teracting spins was also studied to investigate the corre- 
spondence of classical and quantum Liouville dynamics 
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|l3| . But, this correspondence was considered to have 
nothing to do with the decoherence effect. 

The main purpose of this letter is to elucidate that the 
intrinsic decoherence due to the entanglement of subsys- 
tems plays an important role to understand the classical- 
quantum correspondence of composed systems. For this 
purpose, we consider the center of mass motion of the 
delta-kicked rotors composed of two particles of which a 
governing Hamiltonian is given by 

H = ^-(p 2 1 +p 2 2 ) + U(r 1 ,r 2 ) 
Zm 

+fc[cos(r 1 )+cos(r 2 )]^(5(t-iT), (1) 

i 

where r± and r 2 are angle variables with range 27r radians 
and U{t\,t-i) is the interaction potential which confines 
two particles within a distance w. If the confinement 
width w goes to zero, then the system reduces to a usual 
delta-kicked rotor. 

To investigate the center of mass motion, we intro- 
duce canonical transforms, R = {r\ + r<x)j2, r = r\ — r 2l 
P = Pi + P2, and p = (pi — P2)/2. Then, we obtain the 
following Hamiltonian, 



P 2 p 2 

H= + — + U(r) 

2M 2/i W 



Kcos(R) cos(~) **T6(t-iT), 

i 

(2) 



where M = 2m, fi — m/2, and K = 2k. Here, U(r) 
corresponds to the confining potential with impenetrable 
walls at r = ±w. Delta-kicks described by the last term 
of Eq. (0) yield the interaction between the center of 
mass motion and the internal degree of freedom, i.e. the 
motion of the reduced mass p,. 

The time evolution of a wave function r, t) is 

given by simple maps. Between each kick occurring at 
t = iT (i is an integer), the center of mass motion and 
the internal motion evolve independently, so that we ob- 
tain 



ip(R,r, NT + Q- 



exp 
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V(R,r, (N- l)T + 0+). (3) 



The wave functions just before and after the delta-kick 
at t = NT are related by the following mapping: 

^(R,r,NT + 0+) = 



exp 



K cos(R) cos (£) ^{R,r,NT + Q-). (4) 



Combining two maps given in Eq. (Q) and (Q), we nu- 
merically calculate the evolution of the wave function 
*&(R,r,t) for various h and w. For most cases, the 
number of basis states used for the motion of coordi- 
nate R and r are 16384 and 256 respectively, while 



32768 and 512 basis are used respectively for small fi. 
The initial condition is chosen to be the ground state 
i&(R,r,0) = l/^/wcos(irr/2w). 

Classical evolution is obtained from a four dimensional 
map for P, P, r and p derived from Hamiltonian (Q). We 
consider an ensemble of particles which are distributed 
from R = to R — 2ir uniformly with P = 0, and from 
r = — w to r = w with a probability cos 2 (irr j '2w) / 'w with 
p = ±po = ±(wh) / (2w) . From this initial condition, the 
classical evolution is simulated, and the ensemble average 
of the normalized variance of center of mass momentum, 
A 2 = 2 ((P 2 ) - (P) 2 ) /MK 2 is computed. 

As a reference, let us mention the case with a sin- 
gle delta-kicked rotor governed by the Hamiltonian H = 
P 2 /2M + Xcos(P)X;. i '5(^ - iT), which corresponds to 
the case with w = 0. For K = 5, the system is fully 
chaotic, and classical kinetic energy increases diffusively. 
But the diffusion is suppressed quantum mechanically, 
and the classical-quantum correspondence breaks down, 
which is the well-known dynamical localization Q . 

Now, we consider two delta-kicked rotors. The inset in 
Fig. 1 shows the differences of classical momentum vari- 
ances between single and two delta-kicked rotors for var- 
ious w. As w is decreased, the difference between them 
vanishes. Meanwhile, shown in Fig. 1 are the differences 
between the classical and the quantum variances for var- 
ious w with H = 0.07. One clearly sees that the quantum 
localization breaks down when w is increased. The differ- 
ence of A c i and A qm nearly disappears for large w, while 
deviations are observed for small w. The break-down 
of quantum localization, i.e., delocalization, is directly 
visible in the momentum distribution function in Fig. 2, 
where the probability distribution of the center of mass 
momentum P are shown at n = 500 for several h and 
w. The exponential localization is changed into a rather 
broad Gaussian-like profile as we decrease h or increase 
w. Let us note that the delocalization of wavefunctions 
alone is not enough to prove the occurrence of decoher- 
ence. In fact, the delocalization was also observed in 
the study of two interacting particles (TIP) in a random 
potential (IJJl^]. However, it was not attributed to the 
decoherence. 

Next, we consider the reduced density matrix for the 
center of mass motion in order to show that the observed 
classical-quantum correspondence is ascribed to the in- 
trinsic decoherence caused from the interaction with the 
internal degree of freedom. The reduced density matrix 
PR (R 1 ,R 2 ) is given by Tr r (p) = J2 r iK-Ri. r W («2, r). 
As a quantitative measure of decoherence, we calculate 
the linear entropy s; = Tr(pn — p 2 n ) 0. Note that for 
a pure state s; = 0, while for a maximum decoherence 
s; = 1. Figure 3 shows that for large w the entropy si 
rapidly approaches 1, i.e., a maximum decoherence. As 
we decrease w, the entropy si shows rather reduced values 
and slowly increases in time. In fact, an energy level spac- 
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ing of the internal dynamics described by p 2 /2/i + U(r) 
is proportional to 1/w 2 , which means the smaller w, the 
larger the level spacing. For a given perturbation de- 
termined by K it is more difficult to excite the internal 
dynamics in the case of large spacing, which leads to 
effective decoupling between the center of mass motion 
and the internal degree of freedom, and eventually to the 
decrease of decoherence. This is consistent with the pre- 
vious result that the break-down of the localization and 
the classical-quantum correspondence is easily obtained 
for large w, i.e., strong decoherence. 

The cos(r/2) in the last term of Hamiltonian (|^) can 
be regarded as an amplitude noise of the kick onto the 
center of mass motion. As w decreases, the noise is re- 
duced and so is the effect of decoherence. Ott at al. 
studied the effect of noise on the single delta-kicked ro- 
tor, which showed that, for moderate noise and H, the 
diffusion coefficient D qm is proportional to the square of 
the noise level, more precisely D qm ~ v 2 {K/fi) A (y is a 
noise strength) |l6| . For very small fi, it is obtained that 
D qm ~ D c i (D c i is a classical diffusion coefficient). If 
we directly apply this result to our case, we immediately 
obtain the scaling relation D qm ~ (wK/h) 4 for moder- 
ate h since we can approximately regard 1 — cos(r/2) as 
a noise amplitude, and so v ~ w 2 . This scaling is con- 
firmed by numerically calculating D qm as a function of 
{wK/h) in Fig. 4, where % = 0.25, and K = 5. These 
results also confirm the above proposition that the term 
cos(r/2) in the Hamiltonian (||) can be treated as noise 
and thus induces decoherence. 

Finally, we consider the classical and the quantum en- 
tropy production rates of which coincidence has been 
the criteria for the classical-quantum correspondence of 
chaotic systems in studies of decoherence JTo| . Quan- 
tum mechanically, the von Neumann entropy is given by 
S V n — Tr[pn log(pn)]. Due to the difficulty of computing 
the eigenvalues of large matrix, we use a smaller number 
of basis states than that used in the previous simulation, 
4096 for the motion of R, and 256 for r. We choose 
h = 0.07 and can obtain the quantum evolution up to 
7i = 100 with this smaller number of basis. The time 
evolution of the von Neumann entropy consists of two 
different regimes |l7j as shown in Fig. 5. After the initial 
transient, the von Neumann entropy increases logarith- 
mically on time, which exactly corresponds to its classical 
counterpart jll| . Let us consider the phase space (R, P) 
for the center of mass motion. Classical distributions are 
uniform in R and Gaussian in P. The variance of mo- 
mentum increases diffusively. As a result, the classical 
entropy is given by S c i ~ 0.51og(n). In Fig. 5, the quan- 
tum results (dots) show 0.5 log(n) dependence, consistent 
with the classical prediction. 

In summary, we have studied the dynamics of the cen- 
ter of mass motion of two coupled delta-kicked rotors and 
examined its classical-quantum correspondence using the 
momentum distribution and the entropy production rate. 



We show that the correspondence is obtained through in- 
trinsic decoherence which comes from the interaction be- 
tween the center of mass motion and the internal degree 
of freedom. Neither external noise nor an external envi- 
ronment are assumed, here the decoherence results from 
the internal dynamics of the system itself. This decoher- 
ence effect should also be observed in systems composed 
of many particles. Although our study is motivated by 
the dynamics of the classical object composed of many 
constituents, it might also be relevant for molecules com- 
posed of two or many atoms. We hope that similar re- 
sults would be observed in experiments using molecules 
when they are subjected to an external nonlinear force 
depending on the position of the constituent atoms. 



[10] D. Monteoliva and J. P. Paz, Phys. Rev. Lett. 85, 3373 

(2000) . 

[11] S. Nag, A. Lahiri, and G. Ghosh, Phys. Lett. A. 292, 43 

(2001) . 

[12] P. A. Miller and S. Sarkar, Phys. Rev. E 60, 1542 (1999). 
[13] J. Emerson and L. E. Ballentine, Phys. Rev. A 63, 052103 

(2001); J. Emerson and L. E. Ballentine, Phys. Rev. E 

64, 026217 (2001). 
[14] D. L. Shepelyansky, Phys. Rev. Lett. 73, 2607 (1994). 
[15] P. Schmittenckert et al., Phys. Rev. Lett. 81, 2308 

(1998). 

[16] E. Ott, T. M. Antonsen, Jr., and J. D. Hanson, Phys. 
Rev. Lett. 53, 2187 (1984). 

[17] Similar results have been observed in the study of the 
driven double well system coupled to the external heat 
bath In this case, the entropy production rate be- 
comes independent of the diffusion caused by the envi- 
ronment only after some initial transient. 



[1] L. E. Reichl, The transition to Chaos in Conservative 

Classical Systems: Quantum Manifestations (Springer- 

Verlag, New York, 1992) 
[2] W. H. Zurek, Phys. Today 44, No. 10, 36 (1991). 
[3] G. Casati, B. V. Chirikov, F.M. Izrailev, and J. Ford, 

Lecture Notes in Physics Vol. 93 (Springer Verlag, Berlin, 

1979), p. 334. 

[4] G. P. Bermann and G. M. Zaslavsky, Physica (Ams- 
terdam) 91A, 450 (1978); M. Berry et al., Ann. Phys. 
(N.Y.) 122, 26 (1979). 

[5] W. H. Zurek and J. P. Paz, Phys. Rev. Lett. 72, 2508 
(1994); G. Casati and B.V. Chirikov, Phys. Rev. Lett. 
75, 350 (1995); W. H. Zurek and J. P. Paz, Phys. Rev. 
Lett. 75, 351 (1995). 

[6] W. H. Zurek, S. Habib, and J. P. Paz, Phys. Rev. Lett. 
70, 1187 (1993), J. P. Paz and W. H. Zurek, Phys. Rev. 
Lett. 82, 5181 (1999). 

[7] S. Habib, K. Shizume, and W. H. Zurek, Phys. Rev. Lett. 
80, 4361 (1998). 

[8] A. K. Pattanayak, Phys. Rev. Lett. 83, 4526 (1999). 

[9] P. A. Miller and S. Sarkar, Nonlinearity 12, 419 (1999). 



3 




FIG. 1. The differences between classical and quantum 
momentum variances, A ci — A qm , are plotted with M = 1, 
H = 0.25, K = 5, T = 1, and 7i = 0.07. From top to bottom, 
w = 0.0 (i.e. single kicked rotor), 0.1, 0.2, 0.3, 0.4, 0.5, 0.6 and 
0.7. (Inset) The difference between variances of single and two 
coupled delta-kicked rotors, A sing i e — A two for various w are 
shown. From top to bottom, w = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7 
and 0.8. 
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FIG. 2. Momentum distribution functions, (a) h = 0.25, 
w = 0.2 and 1. (b) w = 0.2, h = 0.1 and h = 0.25. Other 
parameters are the same as in Fig. 1. 
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FIG. 4. Quantum diffusion coefficient D qm as a func- 
tion of (wK/h) with h = 0.25. The solid line corresponds to 
Dqm tx (wk/h) 4 . Other parameters are the same as in Fig. 1. 
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FIG. 5. The von Neumann entropies obtained from the 
reduced density matrix p R with Ti — 0.07. For a reference, the 
0.5 log(n) dependence of S c i is represented by the solid line. 
From top to bottom, w = 0.2, 0.4, 0.6, 0.8 and 1.0. 



FIG. 3. Linear entropies for various w. From bottom to 
top, w = 0.1,0.2,0.3,0.4,0.5,0.6 and 0.7. Other parameters 
are the same as in Fig. 1. 
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